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Abstract. The double-mode pulsation of GSC 00144-03031 has been detected when searching for COROT targets. 
A very large dataset composed of 4722 photometric measurements was collected at six observatories in Europe 
and America. There is no hint of the excitation of additional modes (down to 0.6 mmag) and therefore GSC 00144- 
03031 seems to be a pure double-mode pulsator, with a very short fundamental radial mode (P=84 min). From 
uvbyP photometry and evolutionary tracks it appears to be a Pop. I star with M=1.75 Mq, located in the middle 
of the instability strip, close to the Zero- Age Main Sequence. We also discovered other new double-mode pulsators 
in the databases of large-scale projects: OGLE BW2.V142, OGLE BWFV207, ASAS3 094303-1707.3, ASAS3 
000116-6037.0, NSVS 3234596 and NSVS 3324715. An observational Petersen diagram is presented and explained 
by means of new models. A common sequence connecting Pop. I stars from the shortest to the longest periods is 
proposed and the spreads in the period ratios are ascribed to different metallicities (at the shortest periods) and 
to different masses (at the longest ones). 
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1. Introduction 

The double-mode High Amplitude Delta Scuti (HADS) 
stars have been considered for a long time similar to 
double-mode Cepheids, i.e., evolved stars pulsating in the 
first two radial modes only. A remarkable exception is 
constituted by the AC And subclass, composed of stars 
pulsating in the first three radial modes (see Tab. 9 in 
Rodriguez 2003). Very recently, nonradial modes have 
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been found combined with the radial ones, making these 
stars more similar to the low-amplitude, unevolved S Set 
stars which are quite common between the Zero Age 
and the Terminal- Age Main Sequences (Poretti 2003). 
Therefore, HADS stars lost their definition of "simple 
double- or single-mode radial pulsators" . A natural ques- 
tion arises: do all the double-mode HADS stars display 
nonradial modes? To give a first answer to this question we 
investigated in detail the light curve of GSC 00144-03031, 
a new member of the class. Moreover, since other double- 
mode pulsators have been found in large-scale projects, 
it is timely to discuss the entire class by comparing the 
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new, more robust observational results with the theoreti- 
cal models. 

The variability of GSC 00144-03031 was discovered 
when observing the field around HD 43587, a candi- 
date primary target for the asteroseismic space mission 
COROT (COnvection, Rotation and planetary Transits; 
Baglin et al. 2002). The field around this star has been 
explored in the search for secondary targets (Poretti et al. 
2005) and the STARE telescope 1 detected the variability 
of GSC 00144-03031 on the night of 16-17 March 2002 
(JD 2452350). The STARE light curve clearly showed the 
effect of multiperiodicity in this large amplitude variable 
(AV >0.40 mag). 

2. Observations 

Owing to the recent growing interest in high-amplitude 
pulsators, new observations were planned to investigate 
the pulsational content of GSC 00144-03031 in detail. 
The first long subset was obtained at S. Pedro Martir 
Observatory: 249 measurements were collected using the 
simultaneous uvby Danish photometer mounted on the 
1.5-m telescope, on five nights from November 14 to 23, 
2003. The observations were performed with the scope 
of monitoring the brightest COROT targets and there- 
fore the time sampling is not very dense (Fig. ^ upper 
row). HD 43913 was used as comparison star. Later in 
the same winter, dedicated intensive CCD monitoring was 
performed at Athens University Observatory: a 0.40-m 
telescope equipped with a SBIG ST-8 camera and a V 
filter was used on eight nights from January 26, 2004 to 
March 2, 2004. GSC 00144-02970 was used as comparison 
star and GSC 00144-03140 as check star: both are located 
in the same field of view and the monitoring was practi- 
cally continuous (Fig. 2] middle row). A third dataset was 
obtained at Beersel Hills Observatory (Belgium), using a 
0.40-m telescope equipped with an ST10 XME CCD cam- 
era and a V filter. The observations were performed on ten 
nights from February 16, 2004 to March 28, 2004 (Fig. HJ 
lower row). Other less conspicuous subsets were collected 
at SETEC Observatory (0.30-m telescope equipped with 
a ST-8 CCD camera and a V filter, 168 measurements on 
one night), at Sierra Nevada Observatory (OSN; 116 uvby 
data obtained at the 0.90-m telescope on three nights) 
and by MM (0.20-m telescope equipped with a Kodak 
401 CCD camera and a V filter, 483 measurements on 
four nights). The time baseline spans 133 d. 

3. Frequency analysis 

Observations obtained at two different longitudes (Europe 
and America) strongly reduced the daily aliasing effect. 
To perform the various steps of the frequency analysis, we 
used the least-squares power spectrum method (Vanicek 
1971), allowing us to detect one by one the constituents 
of the light curve. After each detection, we refined the 
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Fig. 1. Light curves of GSC 00144-03031 in standard V 
magnitudes. The panels in the upper row show photoelec- 
tric data obtained at S. Pedro Martir Observatory. The 
panels in the middle row show CCD data obtained at 
Athens University Observatory. The panels in the lower 
row are, from left: CCD data obtained at Beersel Hills 
and SETEC Observatories (JD. 2453056); CCD data ob- 
tained at Athens University Observatory and by MM (JD 
2453060); CCD data obtained at Beersel Hills and at 
Athens University Observatories (JD 2453067). 

frequency values by applying the MTRAP code (Carpino 
et al. 1987): the refined values were introduced as known 
constituents in the new search. Such a procedure is partic- 
ularly suitable for a multimode, high-amplitude pulsation 
as it keeps the relationships between the detected terms 
(i.e., 2/i, /i + /2, ■••) locked. Moreover, it does not require 
any data prewhitening as amplitudes and phases of the 
known constituents are recalculated for each new trial fre- 
quency, always subtracting the exact amount of signal for 
any detected term. 

As a first step, we analyzed each subset separately. 
The frequency analysis of the 249 uvby data obtained 
at S. Pedro Martir (SPM subset) clearly detected eight 
terms: /i=17.22 d" 1 , / 2 =22.31 d~\ 2/ X) 3/i, fx + f 2 , 
fi—fxi 2/1+/2, 2/i— /a. The term with the smallest ampli- 
tude is the 2/i —/a one (0.0041 mag). The power spectrum 
does not show any well-defined structure after introducing 
these eight terms as known constituents, even though the 
peaks related to the fi and fa coupling terms (3/i — f 2 , 
3/i + /2,4/i) are embedded among others, without stand- 
ing out clearly. The amplitude of the noise slightly in- 
creases toward low-frequencies (0.0019 mag for / <5 d _1 , 
0.0016 mag for / <10 d" 1 ), while it is quite small at high- 
frequencies (0.0009 mag for 40< / <50 d _1 ). No difference 
was found in the frequency analysis of different colours; 
the rms residuals are 0.0170, 0.0094, 0.0084 and 0.0078 in 
u,v,b and y, respectively: the large scatter in u-light is 
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Table 1. Parameters of the least-squares fits of the measurements on GSC 00144-03031. T =HJD 2453054.1189 
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Fig. 2. Power spectra of the GSC 00144-03031 data. 



due to some instrumental problems with that photomul- 
tiplicr tube. 

The analysis of the 1285 V measurements carried 
out at Beersel Hills Observatory (BHO subset) allowed 
us to detect the same terms as in the SPM subset. 
The amplitude of the noise is 0.0008 mag in the 0< 
/ <100 d _1 interval, with a small increase at low frequen- 
cies (0.0015 mag for / <7d _1 ). The solution with eight 
terms leaves an rms residual of 0.0087 mag. Other combi- 
nation terms are visible in the residual spectrum, slightly 
higher than the noise level. The magnitude differences be- 
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Fig. 3. Light curves of the two independent terms found in 
the time series of GSC 00144-03031. The mean magnitude 
1 / =10.198 has been subtracted from the 4722 datapoints. 



tween the comparison and checks star yield a standard 
deviation of 0.010 mag; the values on individual nights 
range from 0.005 to 0.014 mag. 

The data obtained at Athens University Observatory 
(AUO subset) are the most homogeneous ones. Eleven 
terms have been detected: /i, / 2 , 2/i, 3/i, /1+/2, h~ h, 
2/i +/ 2 , 2/i - f 2 (i.e., those detected in the BH and SPM 
subsets) plus 3/i — /b, 3/1+/2, 4/i (i.e., the terms embed- 
ded in the noise in the BH and SPM subsets). The solution 
with 11 terms leaves an rms residual of 0.0128 mag: it is 
higher than that of previous subsets, but the high num- 
ber of datapoints (2421) allowed us to uncover terms with 
smaller amplitudes, by decreasing the amplitude of the 
white noise (0.0007 mag in the 0-100 d" 1 interval). The 
observed scatter is similar to that observed for the check 
star GSC 0144-03140; the standard deviations range from 
0.010 to 0.018 mag. 
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The three solutions of the V or y light curves are very 
coherent, detecting the same terms with very similar am- 
plitudes and phases (Table ^ a cosine series has been 
used). Therefore, once the least-squares solutions for each 
subset were obtained, the mean magnitudes were aligned 
on the standard value V=10.198. At this stage, the other 
subsets (SETEC, OSN and MM) were added: their sam- 
plings do not allow a reliable frequency analysis by them- 
selves, but combined with the previous ones they supply 
a better frequency resolution. The alignments were made 
by applying the same procedure as above. No systematic 
difference was found between V and y data. Hence, we 
have at our disposal a global dataset composed of 4722 
measurements, which constitutes a solid baseline for the 
careful search for small amplitude modes. 

The analysis of the global dataset evidenced the large 
contribution of the fa term, i.e., 17.22 d~ x (Fig. [21 up- 
per row, left panel). This periodicity is characterized by 
a strong asymmetric light curve since the harmonics 2 fa 
(upper row, right panel), 3/i (middle row, left panel) and 
4/i (middle row, right panel) were detected when going 
further in the analysis. The second contribution in am- 
plitude comes from the /2=22.31 d _1 term (upper row, 
right panel): this term has a sine-shaped light curve, since 
no harmonic was detected. The power spectra also show 
the contribution from the coupling terms between the two 
independent frequencies: the /i + fa and fa — fa terms 
(middle row, left panel), the 2/i — fa and 2 fa + fa terms 
(middle row, right panel) and the 3/i — fa and 3/i + fa 
terms (lower row, left panel). The last power spectrum 
(lower row, right panel) does not show any reliable resid- 
ual peak. The mean noise level is around 0.52 mmag in 
the whole 0-80 d _1 interval, ranging from 0.39 mmag in 
the 50-80 d _1 interval to 0.67 mmag in the 0-10 d _1 one. 

Tabled also summarizes the least-squares solution of 
the global dataset. We also inspected the residual light 
curves to see if there were systematic deviations in some 
part of them, but we found nothing. Formal errors of 
the two independent frequencies were calculated following 
Montgomery & O'Donoghue (1999). Figure |3] shows the 
light curve of each independent term. We stress the fact 
that no additional independent term with an amplitude 
larger than 0.6 mmag was found: therefore, GSC 00144- 
03031 looks like an authentic pure double-mode HADS 
star. 

4. The other double-mode HADS stars 

The class of double-mode HADS stars has been en- 
riched with various new entries provided by large-scale 
projects. In addition to GSC 00144-03031, we discov- 
ered new double-mode pulsators when performing system- 
atic analyses of some public databases: OGLE (Optical 
Gravitational Lensing Experiment; Szymahski 2005 and 
Udalski et al., 1997), NSVS (Northern Sky Variability 
Survey; Wozniak et al. 2004), AS AS (All Sky Automated 
Survey; Pojmanski 2002 and 2003) In some cases, we 
planned new photometric observations. The enlarged sam- 



ple allows us to define better the characteristics of the 
double-mode HADS stars. As shown by the analysis of 
GSC 00144-03031, there are double-mode radial pulsators 
which do not show any hint of nonradial terms, even when 
the amplitude threshold is very small. On the other hand, 
Poretti (2003) showed how the radial double-mode pulsa- 
tion can be also paired to a nonradial one. 

Table El lists the stars where the F and 10 modes have 
been detected. The formal errors on the frequency values 
are usually better than 1-10~ 4 d _1 , which in turn implies 
an error on the ratio values better than 1-10~ 5 . The new 
variables are described below. 

4.1. BW2.V142 and BWLV207 

Moskalik & Poretti (2002) reported on the double-mode 
nature of some stars contained in the OGLE-I database. 
We note that BWl_V109=Macho 119.19574.1169 (Alcock 
et al. 2000). The first least-squares solutions of the 
BW2.V142 and BW1.V207 light curves are given in 
Table Poretti (2000) briefly reported on BW2JVT42 
just showing the light curves on the two terms (note that 
there is a misprint in Fig. 1 about the value of the fa 
term). GSC 00144-03031 and BW2.V142 fill the gap be- 
tween shortest periods (namely those of SX Phe and GSC 
02583-00504) and the group of HADS stars which pulsate 
with periods around 0.10 d. 

4.2. GSC 03949-00811 

The star is catalogued as NSVS 3234596 by Wozniak et al. 
(2004). The frequency analysis of the NSVS data allowed 
us to detect the /i=7.687632 d _1 , / 2 =5.890971 d" 1 , fa + 
fa and fa — fa terms. The amplitude of the 10-mode 
is larger than that of the i^-mode. The scatter is quite 
high (0.034 mag) and it was not possible to evidence other 
coupling or additional terms. The least-squares solution 
of the NSVS light curve is reported in Tab. 01 

4.3. GSC 06047-00749 

Three subsets are available for this star: ASAS 
(ASAS3 094303-1707.3; Pojmanski 2002 and 2003), NSVS 
(NSVS 15731136; Wozniak et al. 2004) and a new one 
obtained by NB (0.20-m telescope equipped with a 
SBIG ST7e CCD and a V filter, 402 measurements on 
four nights in June, 2004). The terms /i = 13.0691 d -1 , 
/2=10.0829 d _1 and fa + fa were detected in all the sub- 
sets. The least-squares solutions supplied a smaller ampli- 
tude for the NSVS subset, since these data are in unfiltered 
light. The solution reported in Table |3| was obtained by 
combining the ASAS and the NB subsets. 

4.4. HD 224852 

The variability of HD 224852=NSV 14800 was discov- 
ered by Strohmeier & Ott (1967) and, independently, by 
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Table 2. Frequencies of the F- and 10 radial modes in double-mode HADS stars 



Star 


F-mode 

r n — 1 T 
[d X ] 


lO mode 
[d J ] 


Ratio 


Reference 


GSC 02583-00504 


19.3347 


25.0066 


0. 


.77318 


Wils et al. (2003) 


SX Phe 


18.19358 


23.37940 


0. 


.77819 


Coates et al. (1980) 


GSC 00144-03031 


17.21752 


22.30760 


0. 


.77182 


This paper 


OGLE BW2.V142 


15.14210 


19.45340 


0. 


.77837 


This paper 


BQ Ind 


12.1951 


15.7686 


0. 


.77338 


Sterken et al. (2003) 


OGLE BW1.V207 


11.68214 


15.09791 


0. 


.77376 


This paper and Moskalik & Poretti (2002) 


AE UMa 


11.62560 


15.03118 


0. 


.77343 


Pocs & Szeidl (2001) 


Macho 116.24384.481 


11.50557 


14.88981 


0. 


.77271 


Alcock et al. (2000) 


Macho 115.22573.263 


10.89872 


14.11001 


0. 


77241 


Alcock et al. (2000) 


RV Ari 


10.73788 


13.89913 


0. 


.77258 


Pocs et al. (2002) 


GSC 06047-00749 


10.08286 


13.06908 


0. 


.77150 


This paper 


Macho 114.19969.980 


9.68320 


12.52957 


0. 


.77283 


Alcock et al. (2000) 


Macho 119.19574.1169 


9.35930 


12.08875 


0. 


.77422 


Alcock et al. (2000); Moskalik & Poretti (2002) 


BP Peg 


9.1291 


11.8329 


0. 


.77150 


Broglia (1959); Figer (1983) 


Macho 162.25343.874 


8.98622 


11.64074 


0. 


.77196 


Alcock et al. (2000) 


AI Vel 


8.96298 


11.59958 


0. 


.77270 


Walraven et al. (1992) 


Macho 128.21542.753 


8.32973 


10.80603 


0. 


.77084 


Alcock et al. (2000) 


Macho 119.19840.890 


7.96395 


10.33167 


0. 


.77083 


Alcock et al. (2000) 


V575 Lyr 


6.86961 


8.96848 


0. 


.76597 


Van Cauteren & Wils (2001; GSC 02118-00297) 


V703 Sco 


6.66838 


8.67921 


0. 


.76832 


Oosterhoff (1966) 


Brh V128 


6.519 


8.497 


0. 


.76721 


Bernhard et al. (2004) 


HD 224852 


6.335592 


8.191974 


0. 


77339 


This paper 


GSC 03949-00811 


5.89097 


7.68763 


0. 


.76629 


This paper 


GSC 04257-00471 


5.753783 


7.514040 


0. 


.76574 


This paper 


VX Hya 


4.47648 


5.78972 





.77318 


Fitch (1966) 



ASAS (ASAS3 000116-6037.0; Pojmanski 2002 and 2003). 
In the framework of our project, the star was observed 
by NB (1671 measurements on nine nights from October 
to December, 2004). The ji=8.1920 d" 1 , / 2 =6.3356 d"\ 
/i + /2, fi — ji terms were detected in the ASAS sub- 
set. The observations carried out by NB are more dense 
and also the 2/2, 2/i + /2 terms were detected. In spite 
of the small differences in the amplitudes between the two 
subsets, the frequency analysis of the combined dataset 
evidenced other coupling terms: 3/i, 2/i — /a, /1 + 2/i. 
The search for new terms was stopped when the noise 
at the lowest frequencies became predominant. The 10- 
mode has an amplitude larger than the F-one and the 
light curves of both periods have a non-sinusoidal shape. 
Table [3] lists the least-squares solution of the combined 
dataset; the mean magnitude has been aligned to that of 
the ASAS subset. The solution with 10 terms leaves an 
rms residual of 0.023 mag; the amplitude of the noise is 
3.3 mmag for / <5 d _1 , 1.6 mmag for 10< / <50 d _1 . 
Also in this case there is no hint of the excitation of ad- 
ditional modes. 

4.5. GSC 04257-00471 

The star is catalogued as NSVS 3324715=NSVS 3361415 
by Wozniak et al. (2004). We could detect the 
/i=5.7538 d"\ / 2 =7.5140 d -1 , 2f x and /1 - f 2 terms 



(Table EJ). The amplitude of the 10-mode is quite small, 
but the NSVS data are in unfiltered light and therefore 
damping is expected with respect to F-light. We note that 
GSC 03949-00811 and GSC 04257-00471 fill an important 
gap toward long periods. 

5. Modelling the double-mode HADS stars 

The period ratios reported in Table are plotted in the 
middle panel of Fig. 01 Most of the stars are grouped 
around log Pq = —1.0 and Pi/Po=0.772. It appears that 
the period ratio decreases toward the longest periods. 
However, some deviating points are observed at both ex- 
trema. In particular, the decrease in the period ratio val- 
ues around log Pq = —0.8 appears much more rapid than 
the one calculated by Petersen & Christensen-Dalsgaard 
(1996; their Fig. 4) using the star sample available at that 
time. Therefore, we calculated a new set of theoretical 
models. 

5.1. The theoretical sequences in the Petersen diagram 

Stellar models were computed with the evolutionary code 
CESAM (Morel 1997). The numerical precision as well as 
the model mesh grid (2000 mesh points approximately) 
were optimized to compute oscillations. The equation of 
state CEFF (Christensen-Dalsgaard & Daeppen 1992) was 
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Table 3. Parameters of the least-squares fits of new double-mode HADS stars. For the stars in the lower part the 
10-mode has an amplitude larger than the F-mode. 



BW2.V142 BW1.V207 GSC 04257-00471 



Term 




Freq. Ampl. 


Phase 


Freq. 


Ampl. 


Phase 


Freq. Ampl. 


Phase 






[d- 1 ] [mag] 


[rad] 


Id" 1 ] 


[mag] 


[rad] 


[d- 1 ] [mag] 


[rad] 


Jl 




15.14210 0.145 


4.84 


11.68214 


0.138 


6.13 


5.75378 0.082 


1.82 






±.00002 ±.006 


±.04 


±.00001 


±.004 


±.03 


±.00008 ±.003 


±.04 


h 




19.45340 0.070 


2.28 


15.09791 


0.022 


0.04 


7.51404 0.033 


1.41 
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used, also including the Coulombian correction to the 
classical EFF (Eggleton et al. 1973). A weak electronic 
screening in these reactions was assumed (Clayton 1968), 
accordingly to the evolutionary stages considered in this 
work. Opacity tables were taken from the OPAL package 
(Iglesias & Rogers 1996), complemented at low tempera- 
tures (T < 10 4 K) by the tables provided by Alexander & 
Ferguson (1994). For the atmosphere reconstruction, the 
Eddington T(r) law (grey approximation) was considered. 

The transformation from heavy element abundances 
with respect to Hydrogen [M/H] into concentration in 
mass Z assumes an enrichment ratio of AY/ AZ = 2 and 
Yp r = 0.235 and Z pr = for Helium and heavy element 
primordial concentrations, respectively. Convective trans- 
port was described by the classical Mixing Length theory, 
with efficiency and core overshooting parameters set to 
uml = ImfRp = 1-8 and d ovel = l ovcv /R p = 0.2, respec- 
tively. The latter parameter is prescribed by Schaller et al. 



(1992) for intermediate mass stars. H p corresponds to the 
local pressure scale-height, while l m and Z ovor represent re- 
spectively the mixing length and the inertial penetration 
distance of convective elements. 

Theoretical oscillation spectra were computed from the 
equilibrium models described above. For this purpose the 
oscillation code Filou (Suarez 2002) was used. The up- 
per panel of Fig. ^ shows various sequences calculated 
from our models. The effects of changing the mass keep- 
ing [Fe/H]=0.00 dex are shown by the three solid lines. 
The other two lines show the effect of decreasing metallic- 
ity keeping M=1.80 M : from [Fe/H]=-0.50 dex (dashed 
line) to [Fe/H]=-1.00 dex (dotted line). 

5.2. The evolutionary tracks 

The physical parameters of HADS stars are sum- 
marized by McNamara (2000). uvby(3 photometry of 
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Fig. 4. Petersen's diagram for double-mode HADS stars. 
The upper panel shows the sequences calculated from our 
models. Solid lines: [Fe/H]=0.00 dex and three different 
masses (from 1.80 to 2.00 M©); dashed line: [Fe/H]=-0.50 
dex and M=1.80 M Q ; dotted line: [Fe/H]=-1.00 dex and 
M=1.80 Mq. The middle panel shows the observed period 
ratios. The lower panel illustrates the matching between 
theory and observations: the solid lines are the theoretical 
sequences calculated for [Fe/H]=0.00 dex, M=1.90 M© 
and [Fe/H]=0.00 dex, M=2.00 M©. The crosses indicate 
low-metallicity (Pop. II) stars; the open symbols (circles, 
squares and triangles) indicate Pop. I stars with different 
masses (1.90, 2.00 and >2.00 Mq, respectively). The filled 
circles indicate the Pop. I stars in which the differences in 
the physical parameters, if any, generate similar period 
ratio values. 



Fig. 5. Evolutionary tracks calculated for [Fe/H]=0.00 
dex, dovcr = 0.2 and four different values for the mass 
(from 1.70 to 2.00 M©). The filled circle indicates the po- 
sition of GSC 00144-03031. The open circles, from bottom 
to top, indicate RV Ari, AI Vel, BP Peg and VX Hya. 



6. Interpreting the Petersen diagram 

The decreasing trends in the theoretical sequences are 
characterized by a long standstill between log Po = —1.30 
and log Po = —0.90 (Fig. upper panel). In this inter- 
val the P\/Pq values are more sensitive to the changes in 
metallicity than to changes in mass. On the other hand, 
small differences in mass produce strong differences in the 
sequences at the longest periods. We can compare these 
predictions with the Petersen diagram as determined from 
the stars listed in Table El 



GSC 00144-03031 (V=10.198, b - y=+0.138, mi=0.194, 
ci=0.849, /3=2.805; SPM and OSN observations described 
above) yields the following parameters: E(,_j / =+0.005, 
M y =+2.00, T off =7590 K, log 5 =3.98, [Fe/H]=-0.05 dex 
(TEMPLOGG method; Rogers 1995, Kupka & Bruntt 
2001). 

The connection between double-mode HADS stars 
is confirmed when considering the evolutionary tracks. 
The tracks have been calculated for [Fe/H]=0.00 dex, 
a value which can be applied to GSC 00144-03031, RV 
Ari ([Fe/H]=0.00 dex), BP Peg (-0.02), AI Vel (-0.15) 
and VX Hya (0.05). Figure [3] shows the positions in 
the log Teff-Mboi plane for these stars. The position of 
GSC 00144-03031 suggests M=1.75 M©; such a position 
is similar to that of RV Ari. Note how the points in Fig. [5] 
describe a narrow instability strip for double-mode HADS 
stars. 



6.1. Stars with P =0.10 d 

The period ratios cluster around 0.772 for log Po = —1.0. 
This point in the Petersen diagram is a well-defined one 
on the basis of the enlarged sample: the theoretical models 
have to match it. We calculated a large variety of models, 
but a satisfactory match is provided only for [Fe/H]=0.00 
dex and 1.70<M<2.00 Mq. These constraints are the 
same obtained from the positions of double-mode HADS 
stars in the log T c ff-Mboi plane (Fig. |5J). We note that 
the changes in metallicity shift the theoretical sequences 
toward higher period ratio values (Fig. 0] upper panel). 
Therefore, the small dispersion of the points related to 
stars with Po=0.10 d suggests a very similar metallic- 
ity (Fig. 01 middle panel). The period ratio value corre- 
sponding to AE UMa (Table EJ) is quite normal (0.7738), 
even though AE UMa is a slightly metal deficient star 
([Fe/H]=— 0.50 dex). However, it is also less massive than 
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other stars (M=1.49 M Q ; McNamara 2000): the two ef- 
fects balance each other and generate the normal period 
ratio value. 

6.2. The shortest periods 

The extreme values observed at the shortest periods re- 
fer to GSC 00144-03031, GSC 02583-00504, SX Phc 
and BW2JV142. If we admit that stars with log P = 
— 1.0 and Pi/Po=0.772 constitute a homogeneous sample, 
GSC 00144-03031 and GSC 02583-00504 can be consid- 
ered as the natural extension toward the shortest period 
of this sample. Indeed, the sequences with [Fe/H]=0.00 
dex and masses in the 1.70-2.00 Mq range match very 
well also these points (Fig. solid lines in the upper 
panel). The high period-ratio values shown by SX Phe 
and BW2_V142 (crosses in the lower panel of Fig. HJ can 
be explained by a lower metallicity (Fig. dotted line in 
the upper panel), which is a well-established fact for SX 
Phe ([Fe/H]=-1.40 dex, McNamara 2000), the prototype 
of Pop. II short-period pulsating stars. The same argu- 
ment holds for BW2_V142, as this star is located in the 
galactic bulge. We can infer that the spread of the period 
ratio values at the shortest periods can be explained by 
different metallicities. 

6.3. The longest periods 

V575 Lyr, V703 Sco, BrhV128, GSC 04257-00471 and 
GSC 03949-00811 are the stars defining the rapid decrease 
at log Pq > —0.90. This decrease is in agreement with the 
theoretical sequences and hence these stars are naturally 
linked to the ones at p)=0.10 d. However, the smaller the 
mass, the more rapid the decrease is (Fig.0J solid lines in 
the upper panel). Therefore, only the M=1.90 Mq models 
are able to match the points related to V575 Lyr, V703 
Sco and BrhV128 (Fig. open circles in the lower panel) 
and only the M=2.0Q Mq models are able to match the 
points related to GSC 04257-00471 and GSC 03949-00811 
(Fig. 01 open squares in the lower panel). 

Two long-period stars show period ratios higher than 
0.770, i.e., HD 224852 and VX Hya (triangles in the lower 
panel of Fig.^J. Such values are admitted by models with 
M> 2.00 Mq only. Indeed, VX Hya is a very luminous 
(see Fig. El and massive (more than 2.30 Mq, McNamara 
2000) star. It is for sure a Pop. I star ([Fe/H]=+0.05 dex) 
and therefore we can rule out a connection with the Pop. II 
stars SX Phe and BW2_V142. 

Therefore the spread in the period ratio values ob- 
served at the longest periods can be explained by dif- 
ferences in the masses. These differences in masses can- 
not be distinguished at log Po < —0.90, where the se- 
quences are almost coinciding. We note that the rapid de- 
crease in the Pi/Pq values for M< 1.80 Mq could explain 
some unusual ratio (as for example Macho 104.20389.1202, 
log Po = -0.80, Pi/P o =0.751; Alcock et al. 2000), with- 



out invoking the excitation of nonradial modes (see Fig. 5 
in Poretti 2003). 

7. Conclusion 

The new observations of GSC 00144-03031 confirm the 
existence of pure radial double-mode pulsators. Thanks 
to large-scale projects it has been possible to discover 
several new double-mode HADS stars and to discuss the 
properties of the class. Detailed models have been cal- 
culated and compared with the observed period ratios 
and the physical parameters determined on the basis of 
uvbyP photometry. Our models suggest that double-mode 
HADS stars constitute a very homogeneous class of vari- 
able stars. The Petersen diagram for these stars can be 
reconstructed using models having [Fe/H]=0.00 dex and 
M=1.70-2.00Mq. Different masses generate similar pe- 
riod ratios for log P < —0.90, but theoretical period ratios 
predict differences for log P > —0.90. More massive stars 
show higher values: this spread is actually observed in the 
Petersen diagram. The sequences with 1.90<M<2.00 Mq 
match most of the observed period ratios in a very sat- 
isfactory way. On the other hand, short-period Pop. II 
stars show slightly higher period ratios owing to their dif- 
ferent metallicity and they clearly deviate from the Pop. I 
sequence. 
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